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Abstract 

We study the critical behavior of a frustrated Blume-Capel (BC) antiferromagnet on a triangular 
lattice by Monte Carlo simulations. For a reduced single-ion anisotropy strength —1.47 < d < we 
find two phase transitions. The low-temperature phase is characterized by the antiferromagnetic 
long-range ordering (LRO) on two sublattices with the third one remaining in a non-magnetic state. 
At higher temperatures there is a critical region of the Berezinskii-Kosterlitz-Thouless (BKT) type 
with a power-law decaying spin-correlation function. For —1.5 < d < —1.47, there is only one 
phase transition from the LRO to the paramagnetic region and the transition is of first order. The 
presence of the BKT phase in the current frustrated BC model is a new feature not observed in its 
non-frustrated counterparts. The values of the decay exponent 77 of the BKT phase corresponding 
to upper and lower temperatures appear to be consistent with the theoretical predictions for the 
six-state clock model. 

PACS numbers: 05.50.+q, 64.60.De, 75.10.Hk, 75.30.Kz, 75.50.Ee, 75.50.Lk 
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I. INTRODUCTION 



It has been shown exactly that in the fully frustrated triangular lattice Ising antiferro- 
magnet (TLIA) with spin 1/2 no long-range order can exist down to zero temperature^ but 
the ground state is critical with the power-law decaying spin-correlation function^. However, 
the situation can change dramatically for larger spin values. The series of studies^"- have 
argued that long-range order (LRO) can occur in the ground state if the spin is larger than 
some critical value. The corresponding spin structure is of the type (1, —1, 0), i.e., with two 
sublattices of opposite magnetizations and one sublattice of zero magnetization. The upper 
bound of this critical value was estimated by the use of Peierls' argument^ as 62 and a more 
precise value was established by Monte Carlo simulations^ as 11/2. Generally, the lack of 
order in frustrated spin systems is due to large ground-state degeneracy and the above stud- 
ies have shown that such degeneracy can be considerably affected by the spin magnitude, 
which can lead to long-range ordering. Nevertheless, the large degeneracy can also be lifted 
by some other perturbations, resulting in long-range ordering even in the highly frustrated 
spin-1/2 system, such as an external magnetic field^""— , selective dilutioniiii^ or inclusion of 
the exchange interactions with further neighborsi^"— . 

It is well known that in the Ising models with spin larger than 1/2 a single-ion anisotropy 
is another parameter that may play a crucial role in their critical properties (see, e.g.,— 
This so called Blume-Capel (BC) model has been intensively studied^i^— mostly on bipar- 
tite lattices, in which case the sign of the exchange interaction is irrelevant to their critical 
properties in the absence of an external field. The model has been confirmed to belong 
to the standard Ising universality clasa^i. However, for an antiferromagnetic BC model on 
non-bipartite lattices we can expect qualitatively different behavior. A frustrated antiferro- 
magnetic spin-1 BC model on a triangular lattice has been investigated by position-space 
renormalization group (PSRG)^'' and transfer matrix^^ methods, and has been found to 
display a finite-temperature antiferromagnetic (AF) LRO of the type (1,-1,0) within a cer- 
tain range of the single-ion anisotropy strength, accompanied with a multicritical behavior. 
Nevertheless, the universality class of the identified second-order phase transition was not 
examined. On the other hand, it is known that a number of frustrated systems violate 
the ordinary universality hypothesis. For example, the spin-1/2 Ising antiferromagnet with 
the frustration arising from the competing nearest-neighbor (NN) and next-nearest-neighbor 
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(NNN) interactions on a square lattice lead to a nonuniversal (or weakly universal) critical 
behavior in which the critical indices of the model depended on the NNN to NN interac- 
tion ratio^"— . Similar behavior was also found in the spin-1 model involving either the 
competing NN and NNN interactions^"^ or positive biquadratic interactions'^-. 

Therefore, the motivation for the present investigations was to study the character of 
the critical behavior of the geometrically frustrated spin-1 antiferromagnet on a triangular 
lattice with the single-ion anisotropy by Monte Carlo simulations. Surprisingly, we found 
that the phase transition from the LRO phase is not second order to a paramagnetic phase, 
as predicted by the PSRG results-^, but of a Berezinskii-Kosterlitz-Thouless (BKT) type to a 
quasi-ordered phase with algebraically decaying spin correlation function, which persists for a 
range of intermediate temperatures between the LRO and paramagnetic phases. This finding 
is unexpected not only because the BKT phase was not found in the earlier investigations22i2i 
but also because no BKT phase was predicted to exist at finite temperatures in the spin-1 
TLIA model^*^. 

II. MODEL AND SIMULATION DETAILS 

We consider the model described by the Hamiltonian 



where Si = ±1, is an Ising spin on the ith lattice site, {i,j) denotes the sum over nearest 
neighbors, J < is an antiferromagnetic exchange interaction parameter, and D is a single- 
ion anisotropy parameter. 

In order to study phase transitions in the present spin system we employ Monte Carlo 
(MC) method. We perform MC simulations on spin systems of the size L^, where L = 
24, 48, 72, 96, and 120. We apply the periodic boundary conditions and the updating follows 
the Metropolis dynamics. To obtain dependencies of various thermodynamic quantities on 
the reduced temperature /cbT/|J|, we use standard MC simulation in which for thermal 
averaging we typically consider up to = 2 x 10® MCS (Monte Carlo sweeps or steps 
per spin) after discarding another A^o = 0.2 x MCS for thermalization. The simulations 
start from high temperatures, using random initial configurations. Then the temperature is 
gradually lowered with the steps /cbAT/| J| = 0.02 (or 0.01 around the critical region) and 
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the simulations start from the final configuration obtained at the previous temperature. In 
order to obtain critical indices, we perform finite-size scaling (FSS) analysis, in which case 
we apply the reweighing techniques^i^ and use N = 10'' MCS. We note that sufficiently long 
simulation times are necessary for the present system, since the integrated autocorrelation 
time at the criticality ranged from r ~ 10^ MCS for L = 24 up to r ~ 10^ MCS for L = 120, 
following the scaling law r oc with the estimated exponent z ^ 2.2. For more reliable 
estimation of statistical errors, we used the F-method^. 

For an antiferromagnet, as an order parameter it is useful to define the staggered mag- 
netization per site as 

= (M,)/L2 = 3( max ( ^ S„ ^ Sj, ^ S^) - min ( ^ S,, J^S^J^ Sk) )/2L', (2) 

i€A jeB keC i€A j^B keC 

where (■ ■ ■ ) denotes the thermal average. Further, the following quantities which are func- 
tions of H or/and Mg are defined: the specific heat per site 



the staggered susceptibility per site 



(3) 



L^ksT ' ^ ' 

the derivatives of the following functions of {Ms) with respect to /3 = l/ksT 

and the Binder parameter (magnetic fourth-order cumulant) 

U = l-^. (7) 

The above quantities are useful for localization of the phase boundaries as well as for de- 
termination of the nature of the phase transition. For example, temperature-dependences of 
a variety of thermodynamic quantities display extrema at the L-dependent pseudo-transition 
temperatures /cbTc(L)/| J|. Thus, for the second-order transition, the critical temperature 
can be estimated from the locations of the peaks of the response functions, such as c and 



4 



Xsi for a given value of L. Then, the observed extrema are known to scale with a lattice size 
as, for example: 

X.(L)(xL^/^ (8) 

Di,(L)ocLV-, (9) 

D2s{L) oc L"\ (10) 

where 7 and v represent the critical exponents of the staggered susceptibility and correlation 
length, respectively. More precise locations of the extrema used in FSS can be obtained by 
reweighing techniques applied to the simulation results performed at the pseudo-critical 
temperature ksT^iL) /\J^^^. 

Furthermore, it is known that in the ground state the sublattice spin-correlation function 
of the TLIA model decays as a power law^: 

{S,S,) (X rT\ (11) 

where r] is the critical exponent of the correlation function. The exponent t] of the model 
with zero single- ion anisotropy has been shown to decrease with the spin value from r] = 1/2 
for spin- 1/2 to zero for spin larger than 11/2, for which the AF LRO occurs^'^. Power- 
law decay of the spin-correlation function is a characteristic of the Berezinskii-Kosterlitz- 
Thouless (BKT) phase^^ and the exponent r] can be estimated by FSS of the order parameter 
which scales as^ 

m,(L) oc L^"^'^. (12) 

Alternatively, it can also be obtained from the staggered susceptibility, which in the BKT 
phase where (M^) vanishes in the infinite lattice size limit is more appropriately defined 
as^^ 

and which scales as 

x'siL) oc L'-\ (14) 

In order to distinguish between the second-order and the BKT transitions, one can employ 
a so called cumulant method^*^, which is based on the behavior of the Binder parameter 
U using the formula 

d{UiL')/dU{L))r^ = {L'/Lf/\ (15) 
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In the case of a second-order transition at the critical temperature Tc the exponent v is finite 
and the correlation length diverges as [(Tc — T)/Tc\^'^ . On the other hand, in the case of a 
BKT transition z/ — > oo, singularities are exponential and the correlation length diverges as 

e = eoexp(a[(r,-T)/rj-^/2). (16) 

Then, if the BKT phase is expected between the long-range ordered (LRO) and the param- 
agnetic (P) phases, the following scaling relations apply: 

m,L'' = /i(L-iexp(ari/2)), (17) 

where h = r]/2, t = (Ti — T)/Ti, T < Ti, and Ti is the LRO-BKT transition temperature, 
and 

x'.L-^ = ML-'exp{at-'/')), (18) 
where c = 2 — ri,t = {T — T2)/T2, T > T2, and T2 is the BKT-P transition temperature. 

III. RESULTS AND DISCUSSION 

Let us first examine the ground state properties for different values of d = D/\J\. For 
d = the ground-state configuration is such that the spins on each elementary triangular 
plaquette sum to ±1. Hence, if we consider a hexagonal plaquette with the antiferromagnetic 
arrangement of the nearest neighbors on the honeycomb backbone, as shown in the inset of 
Fig. [H the central Si spin is "free" and the configurations with any value of Si = ±1,0 are 
energetically equivalent. If c? > the configurations with the values of Si = ±1 are preferred 
and the system behaves like a spin-1/2 Ising model with no long-range orderi. If d < the 
configurations with the values of the central spin S'i = ±1 are suppressed and the ordered 
phase with the antiferromagnetic ordering on the honeycomb backbone and non-magnetic 
states of the central spins, i.e.. Si = 0, can occur. The triangular patterns of the type 
(1, —1, 0) are six-fold degenerate. However, such a state is only favorable for —3/2 < d < 0. 
Below d = —3/2 the energy becomes positive and therefore the non- magnetic state with all 
the spins taking zero value is the ground state. The above cases are summarized in Table [B 

At finite temperatures the system is found to display qualitatively different behavior in 
different regions of the single-ion anisotropy strength. Due to the above arguments, we focus 
on the most interesting region of —3/2 < d < 0. In particular, the behaviors in a broad 
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TABLE I: Types of ground-state spin configurations on triangular plaquettes with the correspond- 
ing energies per site {H) /\J\N , for different anisotropy d intervals. 



d ( 


-00,-3/2) (-3/2,0) 


(0,00) 




(0,0,0) (1,-1,0) 


(1,-1, ±1) 


{H)/\J\N 


-1 - 2(i/3 


-l-d 



region of —1.47 ^ d < Q and a narrow region of —3/2 < d < —1.47 are in more detail 
demonstrated on the selected cases of d = — 1 and d = —1.48, respectively. In spite of the 
PSRG expectation of only one LRO-P phase transition, in the former case, there are two 
anomalies in the temperature variations of various thermodynamic quantities, such as the 
staggered magnetization, the staggered susceptibility and the specific heat, suggesting the 
existence of two phase transitions. From the staggered magnetization dependence in Fig. [1] 
we can observe that as the temperature is lowered some ordering is initiated already above 
^5^2/1^1 ~ 0.5. The phase just below A;b^2/|^| is characterized by a finite value of rris for 
finite L but there is another anomalous increase at fc^Ti/l J| ^ 0.4, which eventually leads 
to the expected AF LRO phase of the type (1, —1, 0) below the temperature fcsTi/l J|. 




FIG. 1: (Color online) Temperature variation of the staggered magnetization m^, for d = —1 and 
different values of L. The inset shows the spin arrangement on a hexagonal plaquette with the 
"free" central spin Si. The up and down arrows denote the spin values +1 and —1, respectively. 
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FIG. 2: (Color online) Log-log plots of the staggered magnetization against the lattice size for 
different temperatures. 
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FIG. 3: (Color online) Temperature variation of the exponent r] and the coefficient of determination 
of the linear fit (measure of goodness of fit) obtained from the FSS analysis, for d = —1. The 
symbols and the error bars respectively represent the mean and the extreme values obtained from 
Eqs. (fT2|) and (fn|). 
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FIG. 4: (Color online) Temperature variation of the staggered susceptibility Xsj for d = — 1 and 
different values of L. 

In the intermediate region between /c^Ti/l J| and kBT2/\J\ the order parameter shows 
only slow decrease with the lattice size L. In Fig. [2] we present the FSS analyses of in the 
temperature region comprising the range {kBTi/\J\, kBT2/\J\). Below kBTi/\J\ ^ 0.4 the 
curves turn upward, which indicates that the values remain finite for L — )■ oo, i.e., the LRO 
phase. On the other hand, above /cbTi/| J| ^ 0.5 the curves turn downward, which indicates 
no ordering in the infinite lattice size limit. Excellent linear fits are obtained within the 
temperatures 0.4 < fc^T/l J| < 0.5, indicating the power-law behavior expressed by Eq. f|T2l) . 
Apparently, the slope and therefore also the value of rj varies with temperature. In Fig. Owe 
plot the value of t], as well as the coefficient of determination as a measure of goodness of 
the linear fit, as functions of temperature using both Eqs. (fT2l) and (fT4l) . At low temperatures 
rris is independent of L and thus r] = 0. Within the interval 0.4 < /c^T/l J| < 0.5 the value of 
r] varies from 0.103±0.001 at fc^T/l J| = 0.4 to 0.219±0.001 at fc^T/l J| = 0.5. We note that 
although the coefficient seems to be constantly equal to one up to almost kBT/\J\ ^ 0.57, 
the inset shows that in fact it starts deteriorating already at /c^T/l J| ^ 0.53. At still higher 
temperatures the linear fit is not appropriate and therefore the finite-size dependences are 
no longer power-law. 

This behavior of the order parameter, along with the observations that the staggered 
susceptibility (see Fig. Hj) diverges for L — )• oo and the peaks of the specific heat (Fig. [5]) 
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FIG. 5: (Color online) Temperature variation of the specific heat c, for d = —1 and different values 
of L. 

are rounded and almost independent of the lattice size above L ^ 48, is very similar to 
that observed in the TLIA model with the ferromagnetic NNN coupling^^, the triangular 
lattice planar rotator model in a six-fold symmetry-breaking field-^-, and the six-state clock 
model^, all displaying a BKT type of the intermediate phase. The existence of the BKT 
phase can be verified by analyzing the behavior of the Binder parameter U according to 
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FIG. 6: (Color online) Binder parameter U{L') plotted against (a) \J{L = 24) and (b) U{L = 72) 
for different values of L' > L. 
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FIG. 7: (Color online) Lattice-size dependence of the Binder parameter U for different tempera- 
tures. Circles, squares and diamonds represent typical behaviors in the LRO, BKT and paramag- 
netic phases, respectively. 

Eq. ( ITSl) . If the intermediate BKT phase exists it should appear in the plot as a line of 



points at which U{L) = U{L'). In Figs. 6(a) and 6(b) the parameter U obtained for L = 24 
and L = 72, respectively, is plotted against those for L' > L. The intersection of the 
resulting curve with the straight line U{L) = U{L') represents a nontrivial fixed point. The 



hump at the values corresponding to the intermediate temperatures observed in Fig. 6(a) 
is apparently a finite-size effect since it disappears when large enough lattice sizes, such as 



those in Fig. 6(b), are considered. Then the data for a rage of intermediate temperatures 



lie on a straight line U{L) = U{L'), implying from Eq. ( |T5|) that u = oo and, therefore, the 
exponential divergence of the correlation length characteristic for the BKT phase. 

The lattice-size independence of the Binder parameter U in the intermediate temperature 
phase is further manifested in the flow diagram of U{L) versus in Fig. ^ li L ^ oo, for 
ksT/lJl < 0.40 the value of U = 2/3 represents the trivial fixed point of the LRO phase, 
for kBT/\J\ > 0.53 the value of U = represents the high-temperature fixed point of the 
disordered phase, and in the intermediate-temperature range for sufficiently large sizes the 
value of U remains constant at a given temperature. 

In order to determine more precisely the values of the LRO-BKT and BKT-P transition 
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FIG. 8: (Color online) Finite-size scaling of (a) the staggered magnetization nis and (b) the stag- 
gered susceptibility, according to the scaling relations (fT7|l and (fTSl) . respectively. 



temperatures, /csTi/l J| and kBT2/\J\, respectively, we further employ scaling relations ((Tj 
and f|T8l) . In particular, we use the values of rj determined by FSS above and tune a common 
value of the parameter a and the transition temperatures A;bTi/|J|, kBT2/\J\ such a way 
that the log- log plots of Eqs. f|T7j) and (|T8l) collapse on the universal curves fi and /2, 



respectively. The plots are shown in Figs. 8(a) and 8(b) The best fit corresponds to the 



values of a = 1.21, fc^Ti/l J| = 0.41 ± 0.01 and fc^Ta/l J| = 0.53 ± 0.01. From the slopes 
of the universal curves /i and /2, —b = —0.06 and c = 1.71, respectively, we extract the 
values of T] = 0.12 ± 0.02 at kBTi/\J\ and r] = 0.29 ± 0.04 at fcsTs/l J|, consistent with 
the values obtained from Eqs. f lT2|) and f|T^ . We would like to point out that these values 
are in a fair agreement with the rj values corresponding to the respective lower and upper 
limits of the BKT transition temperatures in the spin-1/2 TLIA model with competing 
NNN interactions^^, the planar rotator model with six-fold symmetry breaking fields^"^ as 
well as the six-state clock model with both non-frustrated ferromagnetic interactions on a 
square lattice^ and frustrated antiferromagnetic interactions on a triangular lattice^. The 
theoretical prediction for the latter are ri{Ti) = 1/9 and ri{T2) = 1/4. 

The LRO-BKT and BKT-P phase boundaries merge at d ^ —1.47, below which the 
transition changes to the LRO-P type, i.e., between the long-range ordered and the param- 
agnetic phases, and the transition is of first order. The discontinuous character is evident 
from the energy histograms, shown in Fig. [9] for d = —1.48. The histograms are bimodal 
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FIG. 9: (Color online) Energy distributions at the L-dependent pseudo-critical temperatures 
kBTc{L)/\J\ for d = —1.48. Double-peaked structure with deepening barrier between the two 
energy states with the increasing L signals a first-order transition. 

with the dip between the peaks observable already at moderate values of L and approach- 
ing zero as L is increased. Thus our Monte Carlo estimate of the point at which the 
BKT transition lines merge and change to the first-order one [dt, kBTt/\J\) ^ (—1.47,0.35) 
gives the values higher than those for the presumed tricritical point obtained by the PSRG 
method {dt, kBTt/\J\) = (-1.494,0.300)^. Below d = -1.48 the transition temperature 
drops sharply and the tunneling times between the two modes increase enormously. At the 
same time the thermodynamic quantities, such as the staggered magnetization and the 
internal energy e, shown in Fig. [101 start displaying strong hysteretic behavior, associated 
with formation of metastable states, when d is increased and decreased. At sufficiently low 
temperatures the transition point can be approximately located as an intersection point 
of the internal energy dependences in the (i-increasing and (i-decreasing processes. As evi- 

— 1.5, in line with 



denced from Fig. 10(b), for kBT/\J\ = 0.2 the transition occurs at dc ' 
our expectations from the ground-state arguments. 

Finally, in Fig. [TT]we provide a rough estimate of the phase diagram in the d — ksT /\J\ 
plane. The boundaries denoted by the empty symbols are obtained from the specific heat 
maxima, using L = 48, N = 2 x 10^ and three independent MC runs. The low-temperature 
branches denoted by the left- and right-pointing triangles represent the jumps of the energy 
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FIG. 10: (Color online) Variations of (a) the staggered magnetization nis and (b) the internal energy 
e in the decreasing (diamonds) and increasing (circles) single-ion anisotropy d at kBT/\J\ = 0.2. 

(and other quantities) in the (i-increasing and (i-decreasing measurements, respectively, and 
outline the two-phase coexistence region characteristic for first-order transitions below d ~ 
—1.47 (see Fig. [T0|) . Based on the ground-state considerations, the true first-order phase 
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FIG. 11: (Color online) Rough estimate of the phase diagram in the d — /c^T/l J| plane obtained 
from the locations of the specific heat maxima (empty symbols), including a few more precise data 
from the FSS analyses (filled symbols). 
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FIG. 12: (Color online) (a) FSS behavior of the maxima of the direct susceptibility x ^-iid the 
logarithmic derivatives of the first and second moments of the magnetization Di and D2, respec- 
tively, in a log-log plot, and (b) FSS fits of the pseudo-transition temperatures ksT^^^/J, where 
X = x,Di,D2 and v = 0.998, for J > and D/J = -1. 

transition boundary is expected to drop to = — 1.5 at zero temperature. The filled symbols 
represent the transition points determined by the FSS analyses above. It is apparent that the 
locations of the specific heat maxima underestimate the LRO-BKT transition temperature 
but overestimate the BKT-P transition temperature. This behavior is typical also for some 
other systems displaying the intermediate BKT phase^"^"^. 

For the sake of comparison, we also checked the critical behavior of the same model 
but with the ferromagnetic interaction J > and d = —1. In this case we only ob- 
served one anomaly in various thermodynamic quantities, associated with the ferromagnetic- 
paramagnetic phase transition with the power-law scaling of various thermodynamic func- 
tions. The FSS analysis of the direct susceptibility x ^^nd the logarithmic derivatives of the 



first and second moments of the magnetization Di and D2^^i shown in Fig. 12(a), indicate 
that the transition belongs to the standard Ising universality class with the critical expo- 
nents vi = 1 and 7/ = 1.75. Both the critical exponents v = 0.998 ±0.011, 7 = 1.754 ±0.014 
and the transition temperature kBT^/ J = 2.399 ± 0.002, estimated from the scaling relation 
^bT^^xI J = ksTc/ J + aL'^/"^ where ksT^^^j J is the temperature at which the quantity 



X displays a maximum and v is the critical exponent estimated above (see Fig. 12(a)), are 
in a good agreement with the recent high-accuracy MC study results^S. 
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IV. CONCLUSIONS 



We have studied the critical behavior of the BC antiferromagnet on a triangular lattice 
by Monte Carlo simulations and found two kinds of phases within the single-ion anisotropy 
strength —1.47 < (i < 0. Below /c^^i/l^l the system displays the antiferromagnetic LRO on 
two sublattices with the third one remaining in a non-magnetic state. Above kBTi/\J\ for 
a range of temperatures up to kBT2/\J\, the ordering is of the BKT-type with a power-law 
decaying spin-correlation function. For — 1.5 < (i < —1.47, there is only one phase transition 
from the LRO to the paramagnetic region and the transition is of first order. This behavior is 
distinctively different from both the ferromagnetic BC model on a triangular lattice and also 
from a non-frustrated antiferromagnetic model on a bipartite lattice^, which do not exhibit 
the BKT phase and the second-order transition from the LRO to the paramagnetic phase is 
of the standard Ising universality class. However, there are some other systems, such as a 
g-state clock model with g > 4, which have been confirmed to display similar critical behav- 
ior to the present model, featuring the low-temperature LRO, the intermediate-temperature 
BKT and the high-temperature paramagnetic phases. Furthermore, for a selected value of 
the single-ion anisotropy d = —1 the current BC model produced the BKT phase with the 
values of the temperature-dependent exponent t] in the high- and low-temperature limits 
consistent with the theoretical predictions for the planar rotator model with six-fold sym- 
metry breaking fields, f]{T2) = 1/4 and ri{T2) = 1/9^^, as well as those estimated by Monte 
Carlo simulations in the spin-1/2 TLIA model with the ratio of competing NN and NNN in- 
teractions equal to one^^, the planar rotator model with six-fold symmetry breaking fields^ 
as well as the six-state clock model^"^. All these models share the six-fold ground-state 
degeneracy, which we believe is behind the universal behavior at finite temperatures. 

Further, it would be interesting to see how the phase diagram evolves if larger spin values 
are considered. Based on the earlier studies-, for 5* larger than some critical value Sc in the 
ground state the LRO should set-in already at = with the BKT phase transitions still 
occurring at higher temperatures. On the other hand, the increasing spin number is believed 
to change the multicritical behavior in the large \d\ limit of the BC model^i^. Thus, our 
future intention is to extend the present investigations to the systems with larger spin values 
and focus on peculiarities arising from the presence of the geometrical frustration. 
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